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This article summarizes and extends the recent developments in the microscopic modeling 
of the magnetic excitations in cuprate two-leg ladder systems. The microscopic Hamiltonian 
comprises dominant Heisenberg exchange terms plus an additional four-spin interaction 
which is about five times smaller. We give an overview over the relevant energies like the 
one-triplon dispersion, the energies of two-triplon bound states and the positions of multi- 
triplon continua and over relevant spectral properties like spectral weights and spectral 
densities in the parameter regime appropriate for cuprate systems. It is concluded that 
an almost complete understanding of the magnetic excitations in undoped cuprate ladders 
has been obtained as measured by inelastic neutron scattering, inelastic light (Raman) 
, scattering and infrared absorption. 
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1. Introduction 



The origin of high-temperature superconductivity in layered cuprates is still heavily 
\ debated. Nevertheless, it has become evident in the last years that the important 

physics takes place in the two-dimensional copper-oxide planes-^'^. At zero doping, 
these compounds are insulating long-range ordered antiferromagnets which are de- 
scribed by (extended) Heisenberg models^. On doping, they can be viewed as doped 
Mott insulators^. Doping destroys the antiferromagnetic order rapidly. At some finite 
value of doping, the superconducting state is realized. 

The long-range ordered Neel state at zero doping suggests to base a microscopic 
description on magnons which interact with the doped holes. But the holes destroy the 
ordered Neel state giving rise to a spin-liquid state with short range antiferromagnetic 
correlations. Thus, it would be more appropriate to set up a microscopic description 
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in terms of the elementary excitations of the spin-liquid state, which can be fractional 
S = 1/2 spinons^ or integer 5=1 triplons'^. Their interplay with the doped holes 
is expected to be crucial for the understanding of the pairing mechanism. But a 
quantitative theory is so far not in sight. 

There is a different class of systems which are spin liquids already in the undoped 
case due to their lower dimensionality'^'^, namely the cuprate two-leg ladders. Two- 
leg ladders have a finite spin gap^"-^^ and the elementary excitations are triplons 
with 5 = 1. The ground states of these ladders can be viewed as realizations of the 
so-called resonance valence bond (RVB) state, a coherent superposition of singlet 
dimers, proposed by Anderson for the two-dimensional cuprates^'-^^. Thus, the physical 
understanding of the cuprate ladders as doped spin liquids is an important step to the 
understanding of the two-dimensional cuprates. 

There has been decisive progress in the understanding of undoped cuprate ladders 
both in theory^'^'^^'^^'^^'^^"^^ and in experiment^^"^^. It is found that a microscopic 
modeling of the magnetic degrees of freedom has to include a four-spin interaction 
besides the usual nearest-neighbor Heisenberg exchange. The size of the four-spin 
exchange in cuprate ladders is determined to be about 20% of the leading nearest- 
neighbour Heisenberg exchange^^'^^'^°~^^. 

The present article has two objectives. First, it reviews briefly our current under- 
standing of the undoped cuprate ladders. Second, it extends the results for two-triplon 
bound states and spectral densities to the experimentally relevant model including the 
four-spin ring exchange. 

The article is set up as follows. In Sect. 2, we introduce the cuprate ladder com- 
pounds as well as the microscopic model and we give a short overview of the recent 
developments. Then we present a sketch of the method we use to describe the two- 
leg ladder system. In the sequel, the major properties of the two-leg ladder system 
are explored in the parameter regime relevant for the cuprate ladders. In Sect. 4 the 
spectral weights are discussed. Sect. 5 describes the energy properties, namely the 
one-triplon dispersion, the energies of the two-triplon bound states and the positions 
of the multi-triplon continua. The last two parts. Sects. 6 and 7, deal with important 
spectral densities. The focus is laid on the dynamic structure factor which is relevant 
for inelastic neutron scattering (INS). A short summary of inelastic light scattering, 

1. e., Raman spectroscopy, and infrared absorption (IR absorption) is also provided. 
Finally, the article is summarized and concluded in Sect. 8 

2. Cuprate Ladders and the Microscopic Model 
2.1. Materials 

Realizations of two-leg ladders are established in SrCu203 (Srl23) ^^-^^ and 
A14CU24O41 (A14) with A={Sr,Ca,La} ^^■^'^ compounds. The fundamental building 
block in both compounds are edge-sharing copper-oxide plaquettes where adjacent 
copper ions are linked linearly to one another by intercalated oxygen ions. In A14 
corner-sharing plaquettes occur also; they give rise to spin chains. But the correspond- 
ing energy scale is much lower and we do not focus on them here. Sro.4Cai3.6Cu2404i 
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becomes superconducting for high pressure^^ enhancing the interest in this class of 
compounds in particular. 

SrCu203 is the prototype of a system of weakly coupled CU2O3 spin-ladders. The 
relevant orbitals of the copper atoms are the planar Q'x2_y2 orbitals. The superexchange 
runs via the Px or Py orbitals of the intermediate oxygen''^. The interladder coupling is 
significantly smaller because the superexchange via a Cu-O-Cu path with a 90° angle 
is reduced due to the vanishing overlap of the orbitals"^^""^"^. In addition, the interladder 
coupling is frustrated because the spin on one ladder is coupled equally to two adjacent 
spins on the neighboring ladder, which themselves are coupled antiferromagnetically, 
see the topview of a plane of ladders on the right side in Fig. 1. The disadvantage of 
SrCu203 is that it has to be grown under very high pressure so that only small single 
crystals or polycrystals are available''^. This hampers the momentum resolution and 
the counting statistics of scattering experiments, in particular of neutron scattering 
investigations. 

The second class of compounds are the so-called telephone-number compounds 
A14CU24O41 with A={Sr,Ca,La}, see left panel in Fig. 1 Two sorts of planes occur. 
One sort of planes consists of two-leg CU2O3 ladders as described in the last para- 
graph. The other sort consists of corner-sharing CUO2 chains. Both one-dimensional 
structures, ladders and chains, respectively, are oriented along the c-axis. The two 
sorts of planes are illustrated on the right hand side of Fig. 1. It is instructive to 
write A14CU24O41 = (A2Cu203)7(Cu02)io in order to emphasize the existence of 
two structures. This formula assumes a commensurate ratio of seven rungs in the 
ladders matching ten copper sites in the chains. But a closer investigation reveals 
that this is only an approximation; the ratio is in fact incommensurate, see e.g. Ref. 
47 and references therein. 

Ladders and chains form two-dimensional layers which are stacked in b-direction in 
alternation with Ae {Sr, Ca, La } layers, see left panel in Fig. 1. Large single crystals 
of A14 can be grown fairly easily so that inelastic neutron scattering investigations of 
high resolution are possible. 

The distance between two copper atoms is roughly the same in rung and in leg 
direction. The typical order of magnitude of the nearest-neighbor Heisenberg coupling 
is 1000 K. It is similar to the one in the two-dimensional cuprates forming square 
lattices^'^''^^. The coupling between two copper atoms on neighboring ladders is weakly 
ferromagnetic due to the bonding angle of 90°. ^^'^^ The nearest-neighbor exchange 
coupling in the chains is mediated via two symmetric Cu-O-Cu bonds with about 90° 
bonding angle implying a weakly ferromagnetic coupling. Its typical value is an order of 
magnitude smaller than the couplings in the ladders^^. The difference in energy scales 
makes it possible to study the behavior of both substructures, ladders and chains, 
separately. 

Two telephone number compounds, namely Sri4Cu2404i (Srl4) 
and La6Ca8Cu2404i (La6Ca8) are the focus of this work. The system Srl4 is in- 
trinsically doped with 6 holes per unit cell. Using X-ray absorption spectroscopy, it is 
possible to conclude that there are on average 0.8 holes in the ladders and 5.2 holes 
in the chains^^ which can be explained by a higher electronegativity in the chains^^. 
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Fig. 1: (color online) Schematic view of the telephone number compound A14. The 
left panel is a three-dimensional view in the layered material. There are planes of ladder 
structures (plane at the bottom, in the middle and at the top), of chain structure (the 
third and the seventh plane) and the planes containing the A cations. The top views 
of the ladder plane (CU2O3) and of the chain plane (CUO2) are shown in the panels 
on the right side. 

It is possible to substitute Sr^+ by isovalent Ca^+ or by trivalent La-^+ or Y^+. 
The substitution by La^+ reduces the number of holes in the system. The limiting 
case of an undoped sample is reached by the (formal) compound La6Ca8. So far, 
crystal growth has not been able to realize La6Ca8. But there exist crystals up to 
La5,2Ca8.8Cu2404i which can already be viewed as (almost) undoped samples^^ -^®. 

2.2. Minimal Microscopic Model 

We want to focus here on spin isotropic models which do not contain anisotropies. This 
focus is justified by the fact that the anisotropies in cuprate systems are generically 
very small relative to the dominant spin isotropic exchange couplings along linear 
Cu-O-Cu bonds. Most often, Heisenberg exchange couplings are considered in the 
description of the magnetic properties of insulating electronic systems. It has been 
realized early by Dirac^° and it was later introduced in condensed matter physics by 
Thouless^^ that the general expression for the exchange Hamiltonian reads 

00 

^« = EE^«"(-i)''°"na„ , (1) 

n=2 an 

where denotes a permutation operator of the localized spins. The sum runs over 
all permutations labelled by an of n spins. Pa„ is the parity of the permutation 
From the physics point of view, these permutations correspond to the exchange of n 
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spins. The first term represents the well-known two-particle exchange which is usually 
the dominant contribution. 

The importance of the exchange processes with more than two particles was real- 
ized first in the description of the magnetism of solid ^He In -^He, the two-particle 
exchange is small due to a hardcore repulsion of the atoms leading to steric blocking. 
Thus, higher ring exchange processes are important. 

The situation in cuprate systems is different to solid ^He because the dominant 
exchange is the nearest-neighbor Heisenberg exchange. Moreover, the typical geome- 
tries consist of square plaquettes so that the lowest order ring exchange comprises four 
spins. To be precise, we denote the Hamiltonian of the two-leg ladder with additional 
four-spin ring exchange by 

<ijkl> 

Here flyj is the permutation operator of two-particles, r e {1, 2} labels the legs and / 
the rungs. The composite index / = (/', r) combines the leg and the rung information. 
Then the permutation operator of four spins is fljyj^j where < ijkl > restricts the 
sum to four spins on a plaquette. We use the subscript eye from 'cyclic' to indicate 
quantities which refer to ring exchange processes. 

It is convenient to introduce spin operators and to rewrite (2) in terms of a Heisen- 
berg model plus additional four-spin interactions" 

H = J± Si,,S2,/ -|- J\\ St,/Si-,/+1 -|- Hcyc (3) 

/ T,/ 

Hcyc = Jcyc [(Sl,/Si,,+i)(S2,/S2,/+l) -|- (Si,,S2,/)(Si,,+iS2,/+l) 

plaquettes 

— (Si,, 52, (Si,, +152,/)] 

Both Hamiltonians (Eq. 2 and Eq. 3) are the same up to two-spin interactions along 
the diagonal of the spin ladder. It is known that these couplings are small (of the order 
of 0.03J_L so that we restrict the discussion in the following to the Hamiltonian in 
Eq. 3. Note that both Hamiltonians are in use in the literature. The relation between 
the exchange parameters of both Hamiltonians reads 

J± = J1 + \j'cyc (4) 

J|| = J|I + ^Ac (5) 

-Vcyc = Jc'yc ■ (6) 

In agreement with the above general considerations, evidence for significant ring 
exchange processes in high temperature superconductors has been found^"^'^^ very soon 
after their discovery^^. It was observed that the derivation of effective spin models 
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from the underlying Hubbard models implies that the most important corrections to 
the nearest-neighbor Heisenberg exchange are four-spin ring exchange terms. This has 
been substantiated in a sequence of papers considering either the three-band Hubbard 
model^'^'^^'^^ or the single-band Hubbard model^^^^^. The effects of the ring exchange 
on the Raman response have been discussed very early^^'^^'^^. 

For planar cuprates, a first experimental indication for the relevance of the four- 
spin ring exchange was conjectured on the basis of the line shape of IR absorption 
data®^. An unambiguous experimental signature was found by INS in La2Cu04 in 
2001''^. The spin-wave dispersion can be explained only by the inclusion of four- 
spin interactions. The Raman line shape also points towards a significant four-spin 
exchange coupling^*'^^. Recently, it was shown that the magnetic excitations in the 
stripe-ordered phase of Lai.sysBao 125CUO4 are quantitatively described assuming 20% 
ring exchange coupling^^. 

Indications for four-spin interactions in cuprate ladders was found in 1999 by 
Il\l5i7,30,3i La6Ca8. The fits of the experimentally measured one-triplon disper- 
sion by the results for the standard Heisenberg model yielded J\\/J± ~ 2 where Jy is 
the coupling along the legs and J± the one along the rungs of the ladder. The very 
large deviation of this ratio from unity appears astounding in view of the geometrical 
structure of the cuprate ladders which is fairly isotropic. The inclusion of the four-spin 
ring exchange can resolve this discrepancy. But the INS data is not precise enough at 
higher energies to determine the value for the four-spin interaction quantitatively. 

The interaction with light provides a convenient way to obtain information on 
energetically higher lying excitations^^'^^. In particular, bound states can be detected. 
The information on these excitations in turn makes a quantitative determination of 
the ring exchange coupling possible^'*'^''. For La6Ca8 it was found that Jcyc ~ (0.2 - 
Q.25)Jx. J\\/J± ~ 1.2-1.3, and J± = 1QQQ-I200cm~-^ yields a consistent description 
of the available exerimental data^^. Thus the Hamiltonian (2) or the Hamiltonian (3) 
constitute the minimal model for the physics of undoped cuprate ladders. 

The analysis of the shape of Raman lines is another means to determine the 
parameters of the minimal model^^'^^. It reveals that the exchange couplings in Srl4 
and Srl23 differ slightly from those in La6Ca8. Both Srl4 and Srl23 are described 
by X := J||/J_i Ki 1.4 - 1.5 and Xcyc := Jcyc/J± ~ 0.2 - 0.25 The values for Srl4 
are consistent with those for IR absorption data for this compound^°'^^ 

In the preceding paragraphs, we have focussed on the progress in the experimental 
analysis. In parallel, there has been important progress in the theoretical understanding 
of spin ladders. The single-triplon properties were elucidated early^'^'^^. The low- 
lying spectrum in the limit of two coupled chains {J± C J\\) has been analyzed by 
bosonization . Without ring exchange the spin gap A vanishes linearly with Jj_ which 
is confirmed also by quantum Monte Carlo^^. Furthermore, the bosonic field theory 
predicts a singlet state between the triplon state and the continuum. 

The opposite limit {J± > Jy) is the local limit of the system which is also gapped''. 
Combined with the field theoretic and the QMC results we know that no phase 
transition takes place for < x < 00. Moreover, inspection of the interaction of 
two triplons in the limit x -> reveals that they attract each other so that bound 
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States occur^'^"-^'^'''^. The most tightly bound state is always found in the S = 
sector^"* -^^ -^^'^-^. Important progress in high order series expansion for multi-particle 
states^^'^^ '"^"''^ made it possible to follow these bound states to the isotropic regime 
J± PS Besides the energies of the bound states also two-triplon spectral densities 
can now be computed by this approach^°'^^'^^. 

For theoretical completeness, the phase diagram was investigated for all possible 
ratios of the coupling constants. If the four-spin interaction is strongly increased the 
gap in the spin ladder vanishes at a critical value^°. Beyond this value, the character of 
the phase changes. The ground-state is no longer the so-called rung-singlet phase^'^'^-^ 
which is dominated by singlets on the rungs and which is smoothly connected to the 
phase of isolated rungs. The phase for stronger four-spin interaction is characterized 
by a staggered dimerization on the legs of the ladder^°'^^^^®. Apart from these two 
phases, a variety of other phases can be reached for other ratios of Jcyc and J± = Jy 
including negative values^^'^^'^^'^^. We do not discuss the phases other than the rung- 
singlet phase in more detail because the generic experimental systems known so far 
provide spin ladders in the range of the rung-singlet phase only. In view of the ample 
theoretical results for phases beyond the rung-singlet phase it would be exciting if 
experimental realizations for such phases could be found. 



3. Method 

In this part, the essential points concerning the continuous unitary transformations 
(CUTs) are discussed^^'^^^^^ Here we use a particle-conserving perturbative 
(;yj77,78,92 uses the states on isolated rungs as reference. The Hamiltonian of 

a spin ladder including additional four-spin interaction (Eq. 3) is re-expressed by 

^ = HX + XH|| + XcycHcyc (7) 

with X = J\\/J± and Xcyc = -^cyc/Jx as perturbation parameters. 
The CUT is defined by 

9,H(/) = [r?(/),H(/)] (8) 

to transform H{l = 0) = H from its initial form (7) to an effective Hamiltonian 
Hgff := H{l = oc) which conserves the number of elementary triplet excitations on 
the rungs, the so-called triplons'^. Formally, the conservation of triplons is expressed 
by = [H±, Heff]. An appropriate choice of the infinitesimal generator r] is given by 
the matrix elements 

r7„(/)=sgn(c7,-C7,)H,j(/) (9) 

in an eigen basis of H±; the q, are the corresponding eigen values^^. 

The perturbative realiziation of the CUT provides a series expansion for the various 
parts of the effective Hamiltonian: the one-triplon hopping amplitudes and the triplon- 
triplon interaction. We have calculated the hopping amplitudes up to order 11 and the 
triplon-triplon interaction up to order 10 in x and Xcyc. Since the expansion is done 
in the couplings which link the rungs each additional order increments the range of 
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the processes which can be captured. So the order characterizes approximately up to 
which range processes are included in the calculation. The range is measured in units 
of the lattice constant along the legs. 

Next, the observables of interest in the ladder system have to be evaluated in 
order to determine the spectral properties of the two-leg ladder. The observables O 
are transformed by the same CUT 

5,0(/) = [rj(/),0(/)] . (10) 

The four /oca/ operators considered are 

0'(r) = Si,,S2,. = ro' (11) 

0]\r) = S,,,S,,,+i (12) 

= ^ (r-2 + To + + ri'i + Ti") 

(9"'(r) = Sf,-S|, = ri'i+ri"' (13) 
O'^(r) = SI, + Si, = To'^ . (14) 

Here % are local operators changing the triplon number by n. The explicit action 
of these operators is given elsewhere^^. The index / = 1,2 in Eq. 12 denotes the 
leg on which the observable operates. In contrast to the triplon-conserving effective 
Hamiltonian /-/eff, the effective observables also excite or annihilate triplons. The parts 
requiring that an excitation is present before the application of the observable do 
not matter at zero temperature. We have calculated those parts of the effective 
observables which correspond to the creation and annihilation of one and two triplons 
up to order 10. The part creating three triplons is determined up to order 8 and the 
four-triplon parts up to order 7 in x and Xcyc. 

As it is common in physics the symmetry properties of the injected triplons play an 
important role. Here we discuss spin and parity. The first two observables 0' and O" 
excite triplons with total spin zero. Thus they are relevant for optical experiments like 
Raman scattering and IR absorption. The latter two observables inject triplons with 
total spin one. They are important to study the dynamic structure factor relevant for 
INS experiments. 

The parity which is conserved in this ladder system is defined for the reflection V 
about the center-line of the ladder, see Fig. 2. If \n) denotes a state where n rung 
singlets are excited to triplons while all other rungs remain in their singlet state one 
finds V\n) = (-l)"|n), see also caption of Fig. 2. The state \n) is a linear combination 
of many n-triplon states. So no generality is lost in writing 

OM0) = Y^\n). (15) 

n>0 

The parity with respect to V of the observables introduced in Eqs. 11-14 is clear from 
their definition: O'" is odd while O' and O^^ are even. The symmetrized observable 
0" = (Oi + 02)72 is equally even. The parity is conserved in the CUT so that V 
applied to both sides of Eq. (15) requires 

OM-l S"|2n), Oeffeven 
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Fig. 2: The operator V reflects about the depicted axis. A single singlet on a rung 
is odd; a triplet on a rung is even with respect to V. The action of V on the rung- 
singlet ground state is defined to be of even parity V\0) = |0), i.e. we assume an even 
number of rungs. If one singlet is substituted by a triplon in the triplon vacuum |0), 
the resulting state |1) is odd: Vll) = 



An even (odd) parity of Oeff implies that Oeff injects an even (odd) number of triplons 
into the system. The importance of these symmetries will be discussed in more detail 
in Sects. 6 and 7 where the experimentally relevant results are presented. 

The series for the effective Hamiltonian and for the effective observables have to 
be extrapolated in order to obtain reliable results describing the regime of cuprate 
ladders. The extrapolation is performed using the one-triplon spin gap A as internal 
parameter^^'^^'^^. Usually a standard [n, 2] dIog-Pade extrapolant or a [n, 2] Fade 
extrapolant is used for the series expressed in A. The achieved range of reliability 
depends on the quantity under study. The spin gap can be computed down to its 
zero. The matrix elements of the interaction are more difFicult to extrapolate. For 
very low values of Xcyc the ratio x can be taken up to 1.5. For Xcyc around 0.2 the 
value of X should not exceed 1.2 to 1.4. 

4. Spectral Weights 

The analysis of the spectral weights helps to determine the quality of the description in 
terms of a small number of triplons. The chosen description works well, if the weight 
resides mainly in a small number of channels with a small number of triplons . If, 
however, the weight is distributed over a large number of channels the description is 
less well adapted to the problem. The focus is laid here on the spectral weights for 
realistic values of the four-spin interaction. 

4.1. S = 1 

First, the spectral weights for excitations with total spin S = 1 are discussed. The 
local observable considered is the combination of Eq. 13 and Eq. 14 

OS=i(r) = l(oiii(r) + C)iV(r)) . (17) 

This combination corresponds to the spin operator acting on the spin on leg 1. The 
spectral weights were calculated up to order 10 in the one- and two-triplon channel, 
up to order 9 in the three-triplon channel. The total spectral weight /tot is equal to 
1/4 for S = 1 because (S^f = 1/4. 

In Fig. 3a the relative spectral weights /w.rei are shown for Xcyc = (black solid 
lines), Xcyc = 0.1 ((blue) dashed lines) and Xcyc = 0.2 (gray (green) dotted lines) and 
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Fig. 3: (color online) (a) Relative weights //v.rei := //v//tot for the 5=1 operator 
5(0"'(r) + (Eq. 13 and Eq. 14). The In are calculated up to and including 

order 10, 10, and 9 in x for N =1, 2, and 3, respectively. The total intensity /tot = 
Ew=i In is rigorosly equal to 1/4. /sum = {k + h + /3)//tot = kjei + hjei + hjei 
denotes the sum of all plotted contributions; it is to be compared to unity, (b) Relative 
weights for the 5 = operator Si,,Si,,+i (Eq. 12). The 1^ are calculated up to and 
including order 10, 8, and 7 in x for N =2, 3, and 4, respectively. The total intensity 
/tot has been extracted from the 11**^ order result for the ground state energy per 
spin, /sum = (h + h + U)l kot = h.rei + h.rei + U.rei denotes the sum of all plotted 
contributions. 



N 6 {1,2,3}. In addition, the sum of all plotted relative weights Uum = {k + h + 
h)lkot is depicted. The sum is very close to unity. Hence contributions with more 
than three triplons only have negligible weight in the parameter regime displayed. 

The three-triplon sector is extrapolated by a [5, 0] dIog-Pade extrapolant because 
other extrapolants are spoiled by spurious poles. We presume that the depicted ex- 
trapolant of the three-triplon weight is affected by the vicinity of such poles so that 
it overshoots for larger values of x. This conclusion is supported by the observation 
that the sum of all plotted relative weights overshoots in the same fashion as the 
three-triplon weight increases for increasing x. The extrapolations of the one- and 
two-triplon sector are very stable so that it is unlikely that they induce the violation 
of the sum rule for large x. 

The general tendency on increasing the four-spin interactions is the reduction 
of the one-triplon weight and the increase of the weights for two and three triplons. 
Extrapolating the one-triplon spectral weight by dIogPade extrapolants, one recognizes 



I 



I 



February 2, 2008 16:35 WSPC/INSTRUCTION FILE review-ks4-gu3 



Spectral Properties of Magnetic Excitations in Cuprate Two-Leg Ladder Systems 11 

poles which indicate that the one-triplon weight vanishes at the phase transition^^. The 
analysis of h in the vicinity of the phase transition indicates that it stays finite at an 
approximate value of 25% . So the larger part of the weight must have been transferred 
to channels with three and more triplons. We conclude that the description in terms 
of rung triplons becomes inappropriate on approaching the phase transition to the 
staggered dimerized phase^°'^'^~^®. For comparison, we mention that the description 
in terms of triplons on dimers for dimerized and frustrated spin chains appears to be 
much more robust^'^^. 

Fig. 3a shows that for realistic values of the four-spin interaction the one- and the 
two-triplon sector contain by far most of the spectral weight. Therefore, the dynamic 
structure factor for these two channels provides an almost complete description of 
INS experiments on cuprate ladders. The corresponding calculation is presented in 
Sect. 6. 

4.2. S=0 

For the S = case, the local observable O^fp in Eq. 12 is considered. The total 
spectral weight /tot(-f) is calculated reliably from the ground state energy per spin 

eo(x) 

In Fig. 3b the relative spectral weights l^jei are shown for Xcyc = (black solid 
lines), Xcyc = 0.1 ((blue) dashed lines) and Xcyc = 0.2 (gray (green) dotted lines) and 
N e {2,3,4}. In addition, the sum /sum = {h + h + U)lltot of all plotted relative 
weights is depicted. The sum rule is fulfilled quantitatively for all values of Xcyc depicted 
in Fig. 3b. 

Since a triplon has 5=1 there can be no one-triplon contribution /i to the 5 = 
response. The I2 contribution is the leading one. Similar to the 5 = 1 case, the leading 
contribution loses spectral weight upon switching on Xcyc. Simultaneously, triplon sec- 
tors with more than two triplons gain spectral weight. The two-triplon sector contains 
only about 50% of the spectral weight for realistic parameters (xcyc k, 0.2, x > 1). 
This means that there are sizable contributions from channels with more than two 
triplons. This is seen in optical experiments^^ where three-triplon contributions are im- 
portant. In addition, the four-triplon sector matters also, giving rise to finite life-times 
(see discussion below for Fig. 4). 

In summary, the analysis of the spectral weights shows that the description where 
rung triplons are taken as quasi-particles works very well. In the regime of parameters 
which matter for cuprate ladders, the sectors with low numbers of triplons capture 
the essential physics. Sectors with more triplons play a certain, but minor, role. 

5. Energy Properties 

The purpose of this section is to illustrate the energetic properties of the two-leg 
ladder with four-spin interaction in the rung-singlet phase. The influence of the four- 
spin interactions on the relevant energies is investigated. The focus is laid on the 
one-triplon dispersion, on the energies of the two-triplon bound states and on the 
position of the multi-triplon continua. In Fig. 4 the most relevant energies are shown 
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Fig. 4: (color online) Dispersions and multi-triplon energies of the two-leg ladder for 
X = 1. The ring exchange is zero in panels (a) and (b) while it is Xcyc = 0.2 in panels 
(c) and (d). One- and three-triplon energies are shown in the left panels while two- and 
four-triplon energies are depicted in the right panels. Recall that the sectors of even 
and of odd number of triplons do not couple. The gray (cyan) solid lines correspond 
to the lower and to the upper band edge of the multi-triplon continua including two- 
triplon binding effects. The dashed gray (cyan) lines denote the lower band edges of 
these continua neglecting two-triplon binding. 



for X = 1 with Xcyc = and Xcyc = 0.2. The results found are qualitatively generic for 
X = 1.0 - 1.5 and Xcyc = - 0.25. 

5.1. One-Triplon Dispersion 

The one-triplon dispersion can be calculated by very many techniques ranging from 
exact diagonalization over density matrix renormalization and quantum Monte Carlo 
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to diagrammatic techniques and series expansions, see e.g. Ref. 7. For consistency 
we show and discuss the series expansion results here. In the left panel of Fig. 4, the 
one-triplon dispersion is shown for x = 1 with Xcyc = (upper panels) and Xcyc = 0.2 
(lower panels). The one-triplon dispersion diplays a global minimum at /c = vr which 
defines the one-triplon gap A. Additionally, a local minimum exists at /c = 0. This 
minimum can be understood as a crossover property to the physics of spin-chains 
where the dispersion at /c = and at /c = tt are degenerate for the infinite system. 
It is caused by level repulsion between the one-triplon state and the three-triplon 
continuum which approaches the one-triplon state more and more as x — > oc. The 
three-triplon continuum is displayed as gray (cyan) lines in the left panels of Fig. 4. 

The main effect of the four-spin interaction is a global shift of the one-triplon 
dispersion to lower energies. The overall shape of the dispersion is almost unchanged. 
The gap is reduced considerably by about a factor of 2 for Xcyc = 0.2. This property 
in particular explains the importance of the four-spin interaction for a quantitative 
description of cuprate ladders. 



5.2. Two-Triplon Continuum and Bound States 

Regarding two-triplon energies not only the one-particle kinetics of the two triplons but 
also the two-triplon interaction is important^^'^°'^^. Generally, one has a two-triplon 
continuum and two-triplon bound states for a given total momentum. The two-triplon 
energies are shown as black lines in the right panels of Fig. 4. 

There are two regimes. For small momenta k < O.SStt, bound states are absent 
and a wide two-triplon continuum is present. For larger momenta k > O.SStt and 
Xcyc = 0, there are two bound states below the two-triplon continuum. There is one 
bound state with total spin 5 = 1 and one with total spin 5 = 0. The binding energy 
of the 5 = bound state is by a factor 2 to 4 larger than the binding energy of the 
5 = 1 bound state which is expected in an antiferromagnet. The dispersion of both 
bound states has a maximum at /( ss O.Stt and a minimum at k = tt. 

Since a finite four-spin exchange shifts the one-triplon dispersion to lower energies, 
also the two-triplon energies are shifted globally to lower energies. The shape of the 
lower and the upper edge of the two-triplon continuum remains almost unchanged 
by the four-spin exchange because the band edges are solely determined by the one- 
triplon dispersion. 

The effect of the four-spin interaction on the binding energy of the two-triplon 
bound states (black dashed lines in Fig. 4) is most interesting. The four-spin interac- 
tion reduces the attractive interaction between the triplons. Thus the binding energy 
of the bound states is reduced upon increasing Xcyc This effect is present for total 
spin 5 = and for total spin 5 = 1. The effect is more pronounced in the 5 = 1 case 
because the 5 = 1 two-triplon bound state is less tightly bound. As a consequence, 
the 5 = 1 bound state is absent for Xcyc = 0.2. The disappearence of the 5 = 1 bound 
state has interesting effects on the spectral line shape of the two-triplon contribution 
to the dynamic structure factor which will be discussed in Sect. 6. 
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5.3. Multi-Triplon Continua 

Next the relative positions of the multi-triplon continua are investigated. The determi- 
nation of the multi-triplon band edges in combination with the multi-triplon spectral 
weights (see Sect. 4) is important to estimate possible life-time effects which are ne- 
glected in the extrapolation procedure. Recall that the generator (9) orders the states 
according to their number of triplons. This works perfectly if the energy of a state 
with n -I- 2 triplons has a higher energy than all states with n or less triplons''^'^'^. But 
this is not always case, see e.g. Fig. 4. If the correlation between triplon number and 
energy is not valid the state with less triplons can decay into the state with more 
triplons. This hybridization implies an additional broadening of the spectral densities. 
It can be viewed as a finite life-time effect. In this regime, the perturbative CUT based 
on the generator (9) works only approximately. The life-time effects are neglected. 

For simplicity, the n-triplon interactions with n > 2 (Ref. 77) are not included in 
the calculation of the band edges. They are expected to be very small^^. 

The reduction of the one-triplon gap A by the four-spin interaction brings all multi- 
triplon continua closer together in energy. The lower band edge of the three-triplon 
continuum is very close to the one-triplon dispersion for k e [0, 0.47r] (left panels). 
The value for the edge of the three-triplon continuum at /( = is the sum of the 
one-triplon gap and the 5 = two-triplon bound state energy at k = ir. 

The situation for two and four triplons is depicted in the right panels of Fig. 4. 
Both continua overlap strongly for both parameter sets. The overlap of the two- 
and four-triplon continuum is independent of the total spin. This is so because the 
minimum four-triplon energy is always a sum of either two times a scattering state of 
two triplons or the sum of a two-triplon bound state with 5 = and a scattering state 
of two triplons. But the total spectral weight in the four-triplon channel with 5=1 
is smaller than the one with 5 = and the resulting life-time effects are expected to 
be smaller for 5 = 1. Possible effects of the neglected processes are discussed later 
in the sections dealing with spectral densities and with the interpretation of optical 
experiments. 



6. Dynamic Structure Factor 

This part deals with the dynamic structure factor of the two-leg spin ladder with four- 
spin interaction. Results for the one- and the two-triplon contribution which capture 
most of the spectral weight (see Sect. 4) are shown. The obtained spectral densities 
are directly relevant for INS experiments. 

In the discussion of the parity as defined in Sect. 3, we found that the observable 
O'" has an odd parity while the observable O''*^ has an even parity. Thus, the observable 
O'" excites only an odd number of triplons and the observable O'^ excites an even 
number of triplons. It follows that the by far most important one- and two-triplon 
contributions can be measured independently by INS because the two contributions 
have a different parity. In particular, the two-triplon contribution should be accessible 
by experiment because it is the leading contribution with an even parity. 
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Fig. 5: (color online) The one-triplon spectral density li{k) for 0"' is shown for (a) 
Xcyc = and various values of x = {0.25; 0.5; 0.75; 1.0} and for (b) x = 1 and various 
values of Xcyc = {0; 0.1; 0.2}. 



6.1. One-Triplon Contribution 

The most important contribution to the dynamic structure factor is the one-triplon 
contribution. The spectral weight accumulates around k = ir. Upon increasing x, the 
spectral weight decreases for small momenta and concentrates more and more at 
/( = TT as illustrated in Fig. 5a. 

The effect of the four-spin interaction on the one-triplon spectral weight is ex- 
emplified in Fig. 5b for x = 1 with Xcyc = (solid black line), Xcyc = 0.1 (dashed 
(blue) line) and Xcyc = 0.2 (dotted gray (green) line). It can be seen clearly in Fig. 
5b that the effect of Xcyc on the /c-resolved one-triplon spectral weight is similar to 
the effect of x. The spectral weight is reduced at small momenta but it increases 
around /c = tt. These findings suggest that the one-triplon spectral weight close to 
the phase transition to the out-of-phase dimerization is governed by the physics at 
k = TT. Similar results were also found by exact diagonalisation^^. 



6.2. Two-Triplon Contribution 

The two-triplon contribution to the dynamic structure factor is the leading part with 
even parity. Its total spectral weight increases slightly by turning on Xcyc (see Sect. 4). 
The generic relative spectral weight is about 20% to 30% of the leading one-triplon 
contribution for realistic values of cuprate ladders. 

In Fig. 6 the result for the two-triplon contribution to the dynamic structure factor 
is shown for X = 1 and Xcyc = (a), Xcyc = 0.1 (b) and Xcyc = 0.2 (c). As discussed in 
Sect. 5.2 the four-spin ring interaction has a strong influence on the 5=1 two-triplon 
bound state. The attractive interaction is lowered by the four-spin interaction which 
induces interesting changes in the line shape of the two-triplon spectral density. 
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Fig. 6: (color online) The two-triplon spectral density l2{k,u) for O'^ with x = 1.0 
and Xcyc = 0.0 (a), Xcyc = 0.1 (b) and Xcyc = 0.2 (c). The gray (cyan) lines denote the 
lower and the upper edge of the two-triplon continuum. The black lines indicate the 
dispersion of the S = 1 two-triplon bound state. The long-dashed dark lines depict 
the lower edge of the 5 = 1 four-triplon continuum. 
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Fig. 7: (color online) The S = 1 two-triplon bound state energy and the spectral 
weight for x = 1.0 with (a) Xcyc = and with (b) Xcyc = 0.1. Left panels: The 
gray (cyan) line denotes the lower edge of the two-triplon continuum and the black 
line depicts the dispersion of the S = 1 two-triplon bound state. Right panels: The 
momentum-resolved spectral weight of the 5 = 1 two-triplon bound state is shown. 

The dispersion of the 5 = 1 two-triplon bound state is denoted as a solid black 
line in Fig. 6 and in the left panels of Fig. 7. In the right panels of the latter figure 
the /(-resolved spectral weight of the bound state is depicted. On increasing Xcyc, the 
binding energy and the spectral weight of the bound state is reduced. At Xcyc = 0.2, 
the 5 = 1 two-triplon bound state has disappeared because the overall attractive 
interaction has become too small. This finding makes a clear and directly verifiable 
prediction for the results of future INS investigations. While a spin ladder without ring 
exchange shows a bound state also in the 5 = 1 channeP°, the spin ladder with ring 
exchange Xcyc = 0.2 does not. 

The reduction of the binding energy of the two-triplon bound state results in a 
sharpening of the two-triplon spectral density at the lower band edge. The effect is 
most prominent near /c = tt. At the value of Xcyc where the bound state is degenerate 
with the lower edge of the two-triplon continuum, the spectral density displays an 
inverse square root divergence. This corresponds qualitatively to our findings for the 
dimerized and frustrated spin-chain at certain values of the exchange couplings®''. At 
Xcyc = 0.2, where the two-triplon bound state has already dissolved in the two-triplon 
continuum, the very sharp structures near the lower band edge are the remaining 
signature of this divergence. 

The lower band edge of the 5 = 1 four-triplon continuum is depicted by long 
dashed lines in Fig. 6. It can be seen that the overlap of the two-triplon and the four- 
triplon continuum increases with increasing four-spin interaction. Possible life-time 
effects are expected to be of minor importance because the spectral weight of the 
four-triplon continuum increases only very slowly like Uiui) oc (Aw)^^/^ at the lower 
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band edge^'' where Aw measures the distance to the band edge. This extremely slow 
increase results from the hardcore repulsion of the triplons which makes them behave 
like fermions at low energies because the system is one-dimensional. 

7. Optical Experiments 

In this section, two-triplon spectral densities with total spin S = are analyzed. The 
effect of the four-spin interaction on the energy spectrum was already discussed in 
Sect. 5.2. In contrast to the 5 = 1 case studied in the last section, the 5 = two- 
triplon bound state is present also for Xcyc = 0.2 in the range of parameters considered 
here. Thus, the impact of the four-spin interaction on the line shape of the two-triplon 
spectral density is less pronounced than in the 5 = 1 case. 

In the sequel, the complete two-triplon contribution is discussed focusing on x = 1. 
The general trends are the same for other values of x. The second and the third part 
deal with the two-triplon contribution to Raman lines and IR absorption. In all parts, 
results are presented for the observable O' acting on the rungs and for the observable 
O" acting on the legs of the spin ladder. 

7.1. Two-Triplon Contribution 

The two-triplon contribution comprises most of the spectral weight in the 5 = sector 
(see Sect. 4). Nevertheless, the three- and four-triplon spectral weights are sizable for 
X > 1 and especially for finite Xcyc. The latter contributions have to be kept in mind 
when comparing the results to experimental data and in estimating possible neglected 
life-time effects. Recall that the observable O' has only contributions with an even 
number of triplons while the observable 0" includes contributions of odd and even 
number of triplons. 

In Fig. 8 results for x = 1 with Xcyc = 0.0 (a) and Xcyc = 0.2 (b) are shown. The 
upper panels depict the findings for O' and the lower panels depict the results obtained 
for O". Detailed information about the 5 = two-triplon bound state dispersion and 
its /c-resolved spectral weight is presented in Fig. 9. The first effect to note is that 
the overall line shape of the two-triplon spectral density hardly changes in the 5 = 
case when the four-spin interaction is switched on. As in the 5 = 1 sector, there is a 
global shift to lower energies at finite Xcyc resulting from the change of the one-triplon 
dispersion. 

Both observables yield the same line shape for /c = and Xcyc = This 
symmetry holds only at /c = 0; it is broken for finite Xcyc. Nevertheless, even at 
Xcyc = 0.2 the line shape for both observables is very similar for small momenta. A 
more detailed discussion of the k = contribution will be given in the next subsection 
about Raman scattering. For large momenta, however, the responses for the two 
observables differ. 

The next important effect of the four-spin interaction besides the global shift to 
lower energies is the decrease of the binding energy near k = ir. The loss of binding 
energy can be seen clearly in Fig. 9. This effect is accompanied by a reduction of 
spectral weight of the bound state near k = ir, see right panels of Fig. 9. The 
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Fig. 8: (color online) Two-triplon spectral density l2{k,u) with S = at x = 1.0 
with Xcyc = 0.0 (a) and with Xcyc = 0.2 (b). The upper panels correspond to O' and 
the lower panels to O". The gray (cyan) lines denote the lower and the upper edge 
of the two-triplon continuum. The black (blue) lines indicate the dispersion of the 
5 = two-triplon bound state. Long-dashed lines depict the lower edge of the 5 = 
four-triplon continuum. 
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Fig. 9: (color online) Energy and spectral weight of the two-triplon 5 = bound states 
for X = 1.0 with Xcyc = 0.0 (upper panel) and x = 1.0 with Xcyc = 0.2 (lower panel). 
Left panels: The solid gray (cyan) line denotes the lower edge of the two-triplon 
continuum and the solid black line depicts the dispersion of the 5 = two-triplon 
bound state. Right panels: The momentum-resolved spectral weight of the 5 = 
two-triplon bound state is plotted as measured by O' (solid line) and by O" (dashed 
line). 



spectral weight is transferred to the two-triplon continuum which becomes sharper on 
increasing Xcyc. 

The black long-dashed curve in Fig. 8 represents the lower band edge of the 5 = 
four-triplon continuum. The overlap is slightly enhanced by the four-spin interactions. 
Note that the two-triplon bound state is located below the four-triplon continuum for 
all momenta so that no finite life-time effects appear for this part of the two-triplon 
contribution. Furthermore, the two-triplon contribution displays only broad features for 
large momenta and energies above the long dashed curve. Therefore life-time effects 
will have only minor effects on this part of the spectrum because a small additional 
broadening will not change much. 



7.2. Raman Spectroscopy 

This subsection deals with the two-triplon contribution to the non-resonant Raman 
response. The observables for magnetic light scattering in rung-rung (leg-leg) polar- 
ization are the /c = part of O' (O") 22,25,68,69,98,99 |_|g^g Q^|y 3 ^^^^^^ summary of 
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the most important results is given. The reader interested in further details is referred 
to more specialized papers^^'^^'^®. 

Channels with an odd number of triplons are inaccessible by Raman scattering due 
to the invariance of the observables with respect to reflection about the centerline of 
the ladder, see discussion on the parity in Sect. 3. Thus only excitations with an even 
number of triplons matter. The relevant contributions to the Raman response come 
from the two-triplon and the four-triplon sector. The two-triplon contribution is the 
most important one, cf. Sect. 4. 

The line shape does not depend on the observable for Xcyc = because the Hamil- 
tonian is a weighted sum of the two of them H = 0\k = 0) + x(9"(/c = 0) up to 
a global factor^^^s jhus the excited state 0'(/c = 0)|Q) equals -xO"(/c = 0)|Q) 
except for a component proportional to the ground state |0) which does not represent 
an excitation. The Raman response for Xcyc = is discussed in detail elsewhere^^'®^. 
The line shape of 0^{k = 0) and ©"(/f = 0) for finite Xcyc is no longer the same. 
We find that the line shape of 0\k = 0) is almost unchanged by switching on Xcyc. 
In contrast, the line shape of 0"(/c = 0) changes although the effect is not large. 
The overall difference between the line shapes in both polarizations remains small 
for the realistic parameters discussed in this work. The next important effect of the 
four-spin interaction besides the global shift to lower frequencies is a sharpening of 
the dominant two-triplon peak measured by = 0), i.e., in leg-leg polarization. 
In addition, an almost constant plateau is produced for freqencies smaller than the 
two-triplon peak in this polarization. 

7.3. IR absorption 

In this section a short discussion of phonon-assisted IR absorption of magnetic exci- 
tations is presented. The leading infrared-active magnetic absorption is a two-triplon- 
plus-phonon process^^'-^°°'-^°-^. Note that three- and four-triplon processes contribute 
also to the IR absorption signal, depending on the exchange couplings of the lad- 
der, see Fig. 3b. As in the previous subsection, we concentrate here on the most 
important and generic features and we refer to other publications for more detailed 
discussions^"^'^^'^^'^^. 

The two-triplon spectral density l2{k,uj) has to be integrated over all momenta 
weighted by a phonon-specific form factor |fph(/c)p. This yields the two-triplon part of 
the IR absorption /j^^p- The precise form of the phonon form factor |fph(/()P depends 
on the specific phonon involved^°^. Different phonons have different form factors so 
that the sum of several contributions has to be considered. But usually all relevant 
phonon form factors are largest for large momenta at the zone boundary. 

In leg polarization, three features are present in the experimental signals: two peaks 
at lower energy and a broad hump at higher energies. At lower energy, the line shape 
is dominated by the two-triplon 5 = bound state as predicted for small values of 
x.-^^ For X > 0.5, the dispersion of the bound state displays a maximum at /( « 7r/2 
and a minimum at /c = tt, see Sect. 5.2. Both give rise to van Hove singularities in 
the density of states which cause the peaks in /j^^p ■ The spectral weight of the 
bound state has a maximum at k = ir for O", see Fig. 9. Thus both extrema, the 



February 2, 2008 16:35 WSPC/INSTRUCTION FILE review-ks4-gu3 



22 Kai P. Schmidt and Gotz S. Uhrig 

minimum and the maximum, lead to peaks in the two-triplon contribution in leg polar- 
ization which originate from the broadened inverse square root van Hove singularities 
in one dimension. The intensity of the energetically lower peak is larger. The peaks 
are more clearly seen for larger values of x because the difference between maximum 
and minimum of the two-triplon bound state dispersion increases with increasing x. 
The third feature, the broad hump, originates from the two-triplon continuum and 
from three-triplon contributions^'^'^'^. 

In rung polarization, the three features are also present but the intensity is dis- 
tributed differently. The spectral weight of the two-triplon bound state vanishes at 
k = -K (Fig. 9). Thus the energetically lowest van Hove singularity appears only as 
low-energy shoulder. It becomes more pronounced for larger x values. The dominant 
peak results from the van Hove singularity at /c ps 7r/2 where the spectral weight is 
finite. The broad hump induced by the two-triplon continuum is relatively more im- 
portant in this polarization because the two singularities at lower energies carry less 
weight. 

In leg polarization, channels with both odd and even number of triplons con- 
tribute. The analysis of the spectral weights, see Sect. 4, has shown that the three- 
triplon channel contributes significantly. This was to be expected from earlier DMRG 
calculations^"^. Due to the reflection symmetry of the observable, contributions with 
odd parity are absent in rung polarization. Therefore, the next important correction 
is only the four-triplon contribution which carries little weight, see Sect. 4. Moreover, 
its line shape can be expected to be rather featureless and broad, because as many 
as four particles are involved so that momentum conservation does not have a very 
important impact. 

8. Summary and Conclusion 

In this article the spectral properties of magnetic excitations in cuprate two-leg lad- 
der systems have been reviewed and extended to finite ring exchange. A quantitative 
microscopic description comprises two-spin and four-spin interactions. The undoped 
cuprate ladders are always in the rung-singlet phase, i.e. there is an adiabatic connec- 
tion to the state of isolated rung dimers. The elementary excitation is a triplon. The 
exchange couplings generically take the values x = J\\/J± = 12 - 1.5, x^d -li. ~ 0.2 
and Jx = 1000 - 1200cm-i. 

The accuracy of the description in terms of rung triplons in the regime relevant 
for the cuprate ladders was discussed by means of the spectral weights. It has been 
found that the leading one- and two-triplon channels contain the dominant weight. 
A sizable three-triplon contribution is found for total spin 5 = 0. The one- and two- 
triplon channels dominate the physics for realistic values of the coupling constants for 
cuprate ladders. 

Various energies are important: the one-triplon dispersion, the energies of bound 
states consisting of two triplons and the positions of the multi-triplon continua. The 
main effect of the four-spin interaction is a global red shift to lower energies. As 
a consequence, the continua approach each other. The two- and the four-triplon 
continua especially overlap more strongly. The second sizable effect of the four-spin 
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interaction is to lower the attractive two-triplon interaction. Thereby the binding 
energy of the two-triplon bound states is reduced and spectral weight is shifted from 
the bound state to the two-triplon continuum. The 5 = 1 bound state has disappeared 
for realistic values of cuprate ladders due to this effect. 

The dynamic structure factor is dominated by the one-triplon contribution. Most 
of the spectral weight is found at the zone boundaries where the spin gap occurs. The 
spectral weight for small momenta is very low. The two-triplon contribution represents 
about 30% of the total spectral weight integrated over all momenta, frequencies and 
parities. The two-triplon contribution is the leading one with even parity. Thus it 
should be measurable independently from the one-triplon contributions. The shape of 
the two-triplon continuum for total spin S = 1 changes significantly for finite four-spin 
interaction. The disappearance of the two-triplon bound state causes sharp structures 
at the lower band edge of the two-triplon continuum. 

The spectral densities with total spin S = are relevant for Raman scattering and 
IR absorption. The electromagnetic waves modify the exchange couplings on the rungs 
or on the legs of the ladder, respectively, leading to a pronounced dependence on the 
polarization. Raman scattering measures the response at zero total momentum. There 
the spectrum is dominated by the two-triplon peak within the two-triplon continuum. 
The IR absorption displays a structure with three peaks. The two low-energy features 
originate from the two-triplon bound state. The third very broad peak represents 
the continuum contribution. The overall shape of the S = spectral density is not 
changed significantly for realistic values of the four-spin interaction. 

The interplay of recent developments in the theoretical description of gapped 
one-dimensional spin liquids and the improvements in spectroscopic experiments have 
permitted to gain a quantitative understanding of the magnetic excitations and their 
line shapes in cuprate ladder systems. 
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